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0 Filter and method for whitening digitally generated noise. 



0 A digital filter (10) for prewhitening digital noise 
is realized without the need to perform multiplication 
or division operations to weight any of the digital 
samples, instead, the data samples that need weigh- 
ting are bit shifted to achieve multiplication or di- 
vision by a power of two. In the preferred embodi- 



ment, the filtering function is l-z.'8. where z is a 
single sample delay. This function is implemented 
by bit shifting one sample three bits to the right and 
computing the difference between it and an adjacent 
sample. 
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FILTER AND METHOD FOR WHITENING DIGITALLY GENERATED NOISE 



Field of the Invention 



The present invention relates to digital noise 
sources, and more particularly relates to a method 
and apparatus for whitening the noise generated by 
such sources. 



Background and Summary of the invention 

Random noise is used advantageously in a 
variety of applications, such as in exciting electrical 
or mechanical systems under test. In such applica- 
tions, the system's response to the noise stimulus 
can be analyzed to predict the system's response 
to a multitude of other conditions for which individ- 
ual testing is impractical. 

Important in such analyses is the requirement 
that the spectral distribution of energy in the noise 
signal be constant throughout the frequency range 
of interest. If it is not. then subsequent analysis 
cannot determine whether a system's particularly 
strong response at one noise frequency was due to 
a resonance at that frequency or whether the noise 
signal simply had a local maximum at that fre- 
quency component, prompting an exaggerated re- 
sponse. While a noise signal with a uniform spec- 
tral power distribution may be simple to concep- 
tualize, it is often not simple to realize. 

A common technique for generating relative 
low frequency noise (as opposed to microwave 
noise and the like) is to periodically produce a 
digital random number and convert it into analog 
form using a digital-to-analog converter (DAC). The 
analog signal output by the DAC changes at peri- 
odic intervals as new random numbers are pro- 
vided to it. resulting in a time varying "noise" 
signal. This approach is termed the uniformly dis- 
tributed digital random number technique. 

The periodic operation of this random number 
generator affects the frequency spectrum of the 
resulting noise signal. For example, if a new ran- 
dom number is provided to the DAC every four 
microseconds, the resulting noise signal will have a 
nominal bandwidth of (1 4 microseconds), or 250 
kilohertz. Above 250 kilohertz. the power output by 
the DAC begins to roHoff more and more. Even 
below 250 kilohertz. the noise is not perfectly fiat 
but begins to diminish with frequency. These ir- 
regularities are due to a phenomenon commonly 
known as sin{x)/x rolloff, which is present because 
the random numbers are produced periodically and 
the analog samples corresponding thereto are held 
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constant for the periodic interval. {Sin(x)/x is the 
Fourier transform of a constant signal of finite dura- 
tion.) Elimination of this phenomenon would require 
that ideal impulses of voltage be output by the 

5 DAC, a feat that cannot be realized in practice. 

Sin(x).x rolloff is not entirely disadvantageous. 
It tends to limit the power of the noise to a certain 
bandwidth that can be controlled by setting the 
'•sampling frequency" the update rate of the 

w DAC. However, at the corner frequency of the sin- 
(X) X rolloff (250 kilohertz in the above example), 
the response of the noise is already several dB 
down. Ideally, this non-flatness in-band would be 
limited to less than a dB. 

75 This rolloff problem can be overcome by filter- 

ing the noise signal with a filter whose passband 
characteristics compensates for the sin(x)/x droop, 
resulting in an approximately flat response over the 
frequency band of interest. Such filters are termed 

20 "prewhitening" filters and are typically implement- 
ed in digital form for use in the digital portion of the 
noise generator. 

As is known in the art. digital filters operate by 
summing a weighted series of data samples. The 

25 weighting coefficients are selected to achieve a 
desired passband characteristic. Unfortunately, 
digital filters are usually complex and expensive. 
For example, such a filter must usually perform 
several multiplication or division operations per cy- 

30 cle to weight the samples in their proper relation- 
ship. These operations are computationally inten- 
sive and time consuming. Due to the poor perfor- 
mance of conventional microprocessors at these 
tasks, specialized hardware that has been opti- 

35 mized to execute the filtering operations is usually 
required. Such specialized hardware is economical- 
ly impractical to use in many applications, such as 
in simple noise sources. 

The present invention provides the advantages 

40 of digital filtering without the usual high cost and 
complexity. This is accomplished by selecting a 
filter characteristic that can be implemented without 
performing any multiplication or division operations. 
In particular, the present invention implements 

45 a filtering function expressed as 1-z8. where z 
represents a single sample aeiay. The weighting 
coefficients are thus 1 and i 8. The "1" coefficient 
requires no multiplication. Multiplication by the 
"18" coefficient is effected by simply bit shifting 

50 the data sample three bits to the right. The subtrac- 
tion operation is conventional and can be easily 
effected by a simple CPU. The result is a fast and 
simple prewhitening filter that can be easily imple- 
mented without any custom hardware. 

The foregoing and additional features and ad- 
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vantages of the present invention will be nnore 
readily apparent from the following detailed de- 
scription thereof, which proceeds with reference to 
the accompanying drawings. 



Brief Description of the Drawings 



Fig. 1 shows a representative prior art noise 
waveform generated using a uniformly distributed 
digital random number approach. 

Fig. 2 shows the spectral distribution of the 
Fig. 1 waveform. 

Fig. 3 shows the passband response of a 
prewhitening filter according to one embodiment of 
the present invention. 

Fig. 4 shows a prewhitening filter according 
to one embodiment of the present invention repre- 
sented in conventional digital filter topology. 

Fig. 5 shows an apparatus for producing 
prewhitened noise according to the present inven- 
tion. 

Fig. 6 shows a flowchart of a filter process 
according to one embodiment of the present inven- 
tion. 

Fig. 7 shows the spectral distribution of a 
noise source that has been prewhitened using a 
filter according to one embodiment of the present 
invention. 

Detailed Description 



Fig. 1 shows a representative prior art noise 
waveform. As noted, such waveforms are conven- 
tionally generated by periodically computing a digi- 
tal random number and converting the number to 
analog form using a digital-to-ana!og converter. The 
spectral distribution of such a waveform is shown 
in Fig. 2, which clearly illustrates the sin(x)yx func- 
tion. 

As discussed earlier, some prior art techniques 
have sought to mitigate the sin(x)/x effect by pro- 
cessing the digital samples through complex 
prewhitening filters. Other applications use the rip- 
pled noise signal to stimulate the system under 
test, but only analyze the system's response over 
the range across which the noise spectrum is rela- 
tively flat (i.e. the range bracketed and labelled "A" 
in Fig. 2). 

An ideal filter for compensating the sin{x).'X 
noise distribution would exhibit a passband char- 
acteristic of x/sin(x) and would extend the usable 
noise bandwidth of the source up to the first null in 
the noise spectrum. However, this passband cannot 
be realized simply. The present invention uses 
instead a passband function that can be realized 



simply and yet offers first order compensation for 
the sin(x)'X noise function. This function is Fm = Rm 
- Rn-i-'B, where Fn is the filtered noise sample, and 
Rn and Rn-i are sequential random noise samples. 
5 The passband of this function can be determined 
by decomposing the function into its component 
parts. The passband of the origin impulse Rn is 1 
over all frequencies. The passband of the (-1 8) 
impulse at T = 1 is 'C0Sa)-8. which rises from -1 8 to 
10 1 '8 over the frequency range 0 to tt. The summed 
passband is shown in Fig. 3 in solid lines; the two 
decomposed passbands are shown in dashed lines. 
For low frequency noise (i.e. less than about (tt 2)- 
'3. it will be recognized that the gradual increase in 
15 the passband response with frequency closely 
counteracts the gradual decrease in the sin(x):X 
function over the same range. The result is a 
virtually flat noise spectrum over the widened fre- 
quency range labeled 'B.' 
20 Digital filters are conventionally implemented 

by cascading delay stages with output taps at the 
different stages. Delayed signal samples from the 
different taps are applied to respective weighting 
circuits and the weighted samples are summed to 
25 form respective filtered output samples from the 
filter. Fig. 4 shows the present filter function repre- 
sented in this conventional filter topology. As can 
be seen, such an arrangement requires the use of 
at least one multiplier stage to accomplished the 
30 weighting by a factor of -1.8. 

In the present invention, by contrast, this 
weighting operation is accomplished by bit shifting 
the delayed sample 3 bits to the right Bit shifting 
operations are included in the instruction set of 
35 even the most rudimentary microprocessors and 
execute very quickly. The calculation of the dif- 
ference between the delayed, bit shifted sample 
with the current sample is also an operation that 
any microprocessor can perform quickly. 
40 An apparatus 10 for generating and filtering a 

noise sequence according to the present invention 
is shown in Fig. 5. This apparatus is described in 
greater detail in copending U.S. application Serial 
No. 345264 entitled SIGNAL GENERATOR METH- 
45 OD AND APPARATUS, the disclosure of which is 
incorporated herein by reference European Appli- 
cation No . Although not necessary to the 
present invention, illustrated apparatus 10 includes 
a digital mixer for translating the noise spectrum up 
50 from baseband to a desired frequency of interest. 

Apparatus 10 comprises a control microproces- 
sor 12 with associated memory and I/O. a pattern 
memory 14. a pattern menhory controller 16,. a 
digital local oscillator 18, a digital multiplier 20 and 
55 a digital-to-analog converter 22. These elements 
are interconnected by data buses 24 and control 
* lines 26. 

The control microprocessor 12 in the illustrated 
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embodiment is a Motorola 68000 series CPU that 
controls nnost of the functions of a host instrument 
in which the noise source 10 is embodied. Asso- 
ciated with CPU 12 is a memory 28 in which I/O 
and instrument operation routines are stored. Also 
included in this memory are the random number 
generator routine and the prewhitening filter rou- 
tine. Fig. 6 is a flov/chart illustrating the prewhiten- 
ing filter routine. The random number generator 
routine is conventional and many such routines are 
known in the art. 

Pattern memory 14 is a static RAM that can be 
loaded with noise samples produced by CPU 12 in 
accordance with the above-described process. This 
memory has space sufficient to store several such 
noise sequences or may also store other signal 
waveforms, such as a chirp or the like. Provision is 
made for stonng complex pattern data, although 
the noise sequences generated according to the 
present invention are purely real- 
Control microprocessor 12 can be used to con- 
trol all aspects of the pattern memory's operation. 
However, in the preferred embodiment, additional 
memory controller circuitry 16 is provided to re- 
duce the processing burden on the CPU. This 
controller circuitry includes several subcircuits, 
such as a Source Address Register, a DMA Load 
Counter, a RAM Run Counter, a Data Register and 
an End of Cycle Detector. 

The Source Address Register indicates to the 
CPU the base address in pattern memory 14 at 
which to begin loading an incoming noise se- 
quence. The DMA Load Counter indexes from this 
base address to indicate the successive locations 
in the pattern memory in which bytes of an incom- 
ing sequence are to be stored. The Memory Run 
Counter performs the analogous function during the 
reading of the sequence from the pattern memory, 
indexing successive memory locations from which 
sequence samples are to be read. The Data Regis- 
ter holds incoming pattern data that is be loaded 
into the pattern RAM until the RAM is ready to 
allow the transfer. This register permits the RAM to 
output a sequence without delaying the CPU if the 
CPU simultaneously needs to download a new se- 
quence into the RAM. Finally, the End of Cycle 
Detector selects the length of the sequence that is 
to be read from the pattern memory. In the illus- 
trated embodiment, the noise sequence stored in 
RAM can be from 256 to 32768 points long, in 
powers of two. When the End of Cycle Detector 
indicates that a complete sequence has been 
placed back, a variety of operations can be ini- 
tiated. The sequence can start playing back again 
from the beginning, causing a repetitive noise sig- 
nal to be generated. Alternatively, the waveform 
can be turned off. causing a "single shot" noise 
burst. Still further, a signal can be sent to the CPU 



12 to trigger execution of another function. 

Noise sequence data output from the memory 
14 is provided to a digital multiplier 20 that mul- 
tiplies the noise sequence with a data stream cor- 
5 responding to a complex sinusoid from local os- 
cillator 18. Local oscillator 18 is a ROM based 
source that operates in conjunction with CPU 1 2 to 
provide a local oscillator data stream. The mixing 
of these two data streams serves to translate the 
10 baseband noise sequence to a higher frequency. 

The digital data stream output from the multi- 
plier 20 is converted to analog form using a digital- 
to-analog converter (DAC) 22. Since the DAC holds 
each sample for a finite time period (corresponding 
15 to the apparatus's 262.144 Hz sampling rate), a 
sin(x).x term is introduced into the output signal's 
frequency spectrum. This ripple is removed by an 
analog reconstruction filter 34 dedicated to this 
particular function. Filter 34 also serves to suppress 
20 sidebands and harmonics of the desired output 
signal. The filtered analog noise signal is then 
output from the instrument for application to the 
system under test. 

It will be recognized that the process of the 
25 present invention can be used on a continuous 
basts to filter in real time an unending stream of 
digital random numbers as they are produced. In 
one such implementation, the filtering process is 
embedded in the routine that generates the random 
30 numbers itself and is executed as part of the gen- 
eration process. In other embodiments, a limited 
set of random numbers is generated and pro- 
cessed to effect the filtering operation. This latter 
process is illustrated in Fig. 6. A set of 1024 
35 random number samples Rn is calculated and is 
processed to yield a corresponding set of filtered 
samples Frj. 

Having described and illustrated the principles 
of my invention with reference to a preferred em- 
40 bodiment thereof, it will be apparent that the inven- 
tion can be modified in arrangement and detail 
without departing from such pnnciples. For exam- 
ple, while the bit shifting and subtracting operations 
needed to implement the preferred filtering opera- 
45 tion are illustrated as being performed by a micro- 
processor, in alternative embodiments other com- 
ponents may be used. In one such alternative 
embodiment, for example, a shift register is used to 
perform the bit shifting operation. The subtracting 
50 operation is performed by converting the shifted 
sample and the unshifted sample into analog form 
and using an analog operational amplifier to deter- 
mine the difference. The output signal is then al-. 
ready in analog form. Similarly, while the preferred 
55 implementation of the invention uses the filter func- 
tion 1-2.8 (where z is a single sample delay), 
other functions can readily be used v;ith com- 
parable effect. Such other functions include, by 
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way of example, 1/8 - z; 1 - 82: 8 - 1 z: 16 - 2z. etc. 
etc. 

In view of the wide variety of embodiments to 
which the principles of my invention can be ap- 
plied, it should be recognized that the illustrated 
embodiment Is to be considered exemplary only 
and not as limiting the scope of my invention. 
Instead, I claim as my invention all such modifica- 
tions as may come within the scope and spirit of 
the following claims and equivalents thereto. 

Claims 

1. A method of generating a whitened noise 
sequence comprising the steps: 

(a) producing a first random number repre- 
sented in digital form; 

(b) producing a second random number re- 
presented in digital form; 

(c) bit shifting at least one of said samples to 
weight the first and second samples in a one to 
eight ratio; 

(d) computing a difference between the first 
and second random numbers after the aforesaid 
shifting; 

(e) providing said difference as an element 
of the noise sequence; 

(f) producing an additional random number 
represented in digital form; and 

(g) repeating steps (c) through (f) with the 
second random number as the first and with the 
additional random number as the second. 

2. The method of claim 1 in which the repeat- 
ing of steps (c) through (f) is done periodically. 

3. The method of claim 1 in which the bit 
shifting comprises shifting one of said first or sec- 
ond random numbers three places. 

4. A filter (10) for prewhitening a sequence of 
noise signals, comprising: 

input means (12) for receiving first and second 
digitized noise signals; 

bit shifting means (12) coupled to said input means 
for bit shifting at least one of said first and second 
digitized noise signals to weight said signals in a 
one to eight ratio; and 

computing means (12) coupled to the bit shifting 
means for computing a difference between said 
signals after bit shifting and for providing said 
difference as an element of an output, filtered noise 
sequence. 

5. The filter of claim 4 in which the bit shifting 
means includes means (12) for bit shifting one of 
said first or second digitized noise signals three 
places. 

6. The filter of claim 4 in which the input 
means, the bit shifting means and the computing 
means all comprise a microprocessor (12). 



7. In a method of digitally filtering a sequence 
of digital noise samples that includes the step of 
summing a series of weighted data samples, an 
improvement for reducing the computational inten- 
5 sity of the filtering operation, the improvement 
comprising weighting at least one of said digital 
noise samples by bit shifting said sample to 
achieve multiplication or division by a power of two. 
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FIG. 1 
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FIG. 6 
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® Filter and method for whitening digitally generated noise. 



@ A digital filter (10) for prewhitening digital noise 
is realized without the need to perform multiplication 
or division operations to weight any of the digital 
samples. Instead, the data samples that need weigh- 
ting are bit shifted to achieve multiplication or di- 
vision by a power of two. In the preferred embodi- 
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ment, the filtering function is I-2/8, where z is a 
single sample delay. This function is implemented 
by bit shifting one sample three bits to the right and 
computing the difference between it and an adjacent 
sample. 
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